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EXPLOITATION OF DIGITIZED GEOMAGNETIC DATA 


Masahi.sa Suglura 

The lirst mathematical description of the main geomagnetic field 
was given by Gauss^ in 1839. He determined 24 spherical harmonic co- 
efficients based on the data at 84 points read off at 30° intervals of 
longitude along 7 circles of latitude from Sabine's total intensity 
(1837) chart, Barlow's isogones (1833), and Horner's isoclines (1836). 
The field values calctilated from his spherical harmonic expansion were 
then compared with the observations at 91 stations. Since Gauss used 
the 3 components X, \, and Z, the number of values used in his spherical 
harmonic analysis is 252. Following Gauss, Erman and Pedersen' (1874) 
made a main field analysis using ^'alues at 90 points, and thus :»bout 270 
values . 

Adrms^' ^ (1898, 1900) determined 48 coefficients based on the data 
taken at every 10° in longitude and 5° in latitude from Sabine's and the 
Admiralty charts. His data set thus consists of about 3800 values. The 
subsequent main field analyses made by Schmidt'''’'' ’’’ (1889, 1895, 1898), 
Dyson and Furner'' (1923), Bauer® (1923), and other workers used data 
sets in which the number of values is of the order of a few thousand. 
Vestlne's main field analysis^ (1947) used data of a similar size foi. 
the spherical harmonic analysis, but the preparation of the data was 
more elaborate than In the previous attempts. 

The development of high speed computing machines brought a sudden 
change in the scope of the main field analysis. Cain et al.^^ ( 1967) 
derived 120 spherical harmonic cotfflcients and their first and second 
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time derivatives from a sample of all magnetic survey data available for 
the interval from 1900 to 1964 and the total field data obtained by the 
OGO 2 satellite. Their data set contained approximately 77,400 obser- 
vations of one or more components each; the total number of values used 
was roughly 150,000. The data for 1965 included 22,300 OGO 2 observations. 

The scope of the data obtained by satellites is enormously large. 

For Instance, the number of observations provided by the OGO 2 and 4 
satellites is as follows: 


Satel lite 

Year 


Approximate No. 
of observations 

OGO 2 

1965 


2,500,000 


1966 


9,000,000 


1967 


800,000 

OGO 4 

July-Dee . 

1967 

14,500,000 


Jan. -May 

1968 

12,000,000 


Now we turn to Sq analyses. For his classical analysis of Sq, 
Chr^pman^** (1919) used 21 observatories for the years 1902 and 19C5. The 
number of values used is thus 3 (components) x 21 *ob.servator ies) x 
24 (hours) x 365 (days) - 551,880 per year; and hence the total number 
is about 1,100,000 for the two years. 

In Vestine's analysis of the Second Polar Year data he used 55 
stations, making the total number of values approximately 1,400,000. 
Tliese analyses were made manually with desk calculators. 

The analysis of Sq has also seen a drastic change in scope with 
the advent of electronic computers. Matsushita and Maeda^*”^ (1965), 
using the data obtained at 69 stations on 45 quiet days during the IGY, 
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derived Sq current system.s for three longitude sectors and for thri.e 
seasons. The number of data values was of the order of 200,000; but 
the number of current systems presented was 40, which represents an in- 
crease by a factor of 10 compared with previous papers. To obtain these 
results they performed spherical harmonic analysis 40 times, 

Sugiura and Hagan' ^ ^ ( 1967) constructed a motion picture to de- 
scribe a continuous variation of Sq by a combined use of a computer and 
a cathode ray type plotter. They used the IGY data from 66 observatories 
for two magnetically quiet periods of five days each and performed a 
spherical harmonic analysis for each of the 240 hours to represent in- 
stantaneous Sq . Then Sq potential contours, both for the external and 
the internal parts, were drawn by machine. To make a motion picture, a 
smooth transition from one contour map to the one for the next hour was 
achieved by Interpolating contours. In the motion picture the Sq potential 
contour map is given at 3 minute intervals. Thus the total number of 
contour maps is 4800. 

Price and Stone^^' in England are developing a new technique to analyze 
Sq automatically with a computing machine on a routine basis. 

In the efforts of Sq analysis by Matsushita and Maeda'’’ , Sugiura 
and Hagan'^'^^ and Price and Stone',® tabulated hourly values had to be 
punched on cards to be fed into the computer. This is by no means a 
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trivial, task. 
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For several years the U.S. Coast and Geodetic Survey and the Goddard 
Space Flight Center have jointly made great effort to digitize magnet* - 
grams^’' . The project was in part financed by the National Science Foun- 
dation. In this program, magnetograms are scaled by a machine operated 
manually, and the scaled results, after being converted to values in 
proper units, are stored on magnetic tape. Using such digitized data, 
Heppner^^ (1969) has made a motion picture depicting changes in the hori- 
zontal magnetic vectors at high latitude observatories. 

Usefulness of geomagnetic activitv indices is now well recognized. 
Thus, for instance, the Dst inde.x, representing mainly the ring current 
intensity, has been published since IGY e.xcept for a gap of 1959 and 
1960^ . Until the digitization program began in 1964 the hourly 
values of H had to be punched on cards, but with the 2.5 minute scalings 
available on magne tic tape no punching task is needed. 

The derivation of the Dst index, simple as it may appear, does in- 
volve an appreciable amount of work. First, for each station contributing 

0 

to the index the secular change has to be evaluated. This is done in 
the present scheme by taking annual means for quiet days and by expressing 
the secular change in power series in time, thus enabling the 'base line' 
to be determined at any given hour. Using annual means of ten quiet 
days the change in H at Hermanns amounted to 670v *n 9 years from 1954.5 
to 1963.5! Though this represents an exceptionally large variation, 
accurate evaluations of the secular changes are essential for an accurate 
determination of Dst. After the secular change is removed the Sq variation 
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must i)t' determined and likewise removed from the dfU.i. This is a dU’fi- 
cult task, because on a disturbed day Sq cannot bo determined without 
knowing the disturbance variation which we want to derive. At present 
the following method is being used; for each station Sq for each yi’ar 
is represented by a double Fo»irier series with local time and a month 
number (in decimals) - expressing seasc>nal variation - as cwo variables 
so that Sq can be evaluated at any local time on any day of the year. 

The Fourier coefficients are determined by monthly mean Sq as obtained 
from five International quiet days. 

All these calculations are now programmed for .paehine comput icion . 

The approximate number of 2.5 minute scalings used for the Dst derivation 
is 210,000 per year per station. Without a higli speed computing mi>chin>.' 
and a high speed plotter it is virtually impossible to produce the Dst 
index on a routine basis. This applies to the auroral electro jet index, 

.\E . If we require 10 stations for the deriv. :ion of AE , approximately 
i million data values per year have to be dealt with. 

The size of the data required for geomagnetic studies is increasing 
tremendously. At the rate of one reading every 2.5 minutes the number 
of values for the three components is about 6JO.OOO per year per station. 
For 50 stations the number becones 350 million for one solar cycle. 

Even with hourly values the corresponding number is 14.5 million. 

The cost of digitization for Dst and AE is approximatt* ly $100,000 
a year, or about $6,000 a year per station on the average (ignoring dif- 
ferences in cost because of the differences in the complexity of variations 
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at Uw and high latitude observatories). in 5 years this amounts to 
$30,000 per station. An estimate for the production in U.S. of a set 
of fluxgat;* magnetometer, a tape recorder, and other necessary aqulpiuenc 
wliich will produce digital daca automatically is on the order of $30,000 
including labor. Once the magnetometer-recorder system is Installed, 
the sum of $6,000 a year, which would be the average digitization cost 
had the conventional magnetometer ueen operated, would seem to be adequate 
for maintenance of the equlnment. There would of course bt. a saving of 
cost in operating an automatic observatory because of the reduced require- 
ment in manpov/er in the dally care of the record and the preliminary 
scaling tliat is often done at an observatory. 

In addition there are hidden costs in analyzing convent lonrl magneto- 
grams by scientists. These include scalings (often duplicated), reproducing 
rotord.s on different scales, punching the scaled data, labor spent by 
scientists tlu'mselves and accompanying administrative expenses. 

Ia'L UK .summarize at this stage the above overall discussion of the 
uigitii-'Lion of the conventional magnetograms and the automatic iiu> t ic 
observatory. First, there is a growing need for digitized magnetic data 
in a machine-readable form, and tliis need is rapidly becoming a necessity. 
Secondly there may actually be a case of a cost reduction in the long 
run in the observation and analysis of certain aspect, if not all, of 
geomagnetic variations if the conventional observatory is modernized by 
Installing a magnetometer witii an automatic recording system with digitized 


output . 
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On tlK.* seci'nd jioLnt above, I w'ould like Lo 'irgc the lACA to inaKC a 
careful study of the financial questions us well as on the instrumental 

problems. It Is a responsibility of Llie scientists to tlu* soc'ety to 

< 

assure that the funds are being spent in the most tlficlent way. Scienti- 
fic merits and priorities must be carefully assessed r.'garctless of 
wliether the conventional magnetic observations are continued or a new 
system is to be introduced. 

I will now make a little more teciinical comparison of the convent iona 1 
and automatically digitized geomagnetic records, for any extensive 
digital data analysis the advantage of automatically digitized records 
is obvious because of (a) the capability of direct input into a computing 
machine and (b) Llie readiness tor display of tl 2 data in any desired fonn 
h\ a plotting machine of cither .nechanical or cathode ray type. For 
synoptic studies also, automatically digitized records have advantages; 
for instance, (a) for tlie drawing (jf disturba :e field vectors or equi- 
valent current vectors the conventional magnetograms would n-quire a 
scaling; (h) Cor a compari.iOn of records from many stations the cenvenL ic-ua I 
magnetograms often have to be scaled and replotted because of different 
magnetometer sensitivities and different scales both in time and in 
ordinates. inconveniences with the conventional magnetograms in all 
these matters have hindered a rapid progress. 

There are, however, arguments against an automatic magnetometer 
system. But tliese arguments often have pitfalls and fallacies. For in- 
stance, it is sometimes said that the traces in the photographic magneto- 
grams show how rapid the variations, for example, in sudden commencements 
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(SC's), arc, and tUaL tiiis quality will be lost in records ui^itL>u-d at 
frequently su^^ested sanipling intervals. Although what is said luTe is 
true, this ari;unienL cannot claim an advantage of the ‘convent iona i 'Tiau- 
netograms on a scientific ground. The rv*asons for this are as fellows. 

First, the normal magneto grams are not intended and not suited for 
studying variations of time scale ■' 1 minute. One must guard himself 
against an Illusion cau.sed by the pleasing appearance of beautifully 
processed photographic magnetograms. This aesthetic enchantment may 
not lead to anything significant scientifically. For instance, studies 
have been made on possible differences in the arrival tines of SC's at 
different local times. However, no definite conclusions, to my knov.'- 
ledge , have been drawn on tlie question. Considering the complex structure 
of the plasma and tlic mrf:,netic field in the magnetosphere, it seems un- 
likely that tlie propagation of the SC disturbance can be studjed in such 
a detail that the thinness of the magnetogram traces becomes a signifi- 
cant factor. The OGO i and 3 satellite observations witli a great time 
resolution have indicated that the SC-associated magnetic variations are 
indeed very complex in the magnetosphere. What is observed on the ground 
is an integrated effect of the magnetospher ic response to a discontinuity 
in the solar wind. In addition to the principal hydromagne tic effects 
there must be variations arising from currents along the field lines 
and in the ionosphere due to the deformation of the magnetosphere; and 
all these variations are further modified appreciably by tne induced 
currents in the ocean and under the ground. 


TIu'rr may arise a need fur cunvc‘nt it)iij i maj^.!u torrams fo lu'tvu "a 
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little lie tter time resolution" than in sampled il i t i/ed data. To li.- 
prepared for such a need there could he a >;roiip of eonvi tit ion.i ! oli-;i’i - 
vatories maintained as in the* past. T.n some areas ot t!ie world it r.iplit 
not be i;asy to convert tlie convenlional t.«bsei vat»n res into auton.itit. 
observatories for various reasons. Nor would it not be advisable to 
convert every existing observatory into one equipped with an autom.'^tic 
system. There could be two groups of observatories conventional .niid 
automatic serving for different purposes. 

However, the cpiality tif the magnctog.rnms that is consii'f red to lie 
advaritagcous aliovc ccrfisti Lu L e « also a liability. During sc'vere distu.r- 
baiices v/e so often lose track of the traces bec^u.ie .of Lbi'ir tliinness 
in addition to tlu-ir crossing each other repeatedly. Imleinl this is one 
of the strongest argumi-nts against the conventional magne Lo>;ra'ns . 

Atii ther argument frequently presented against an automatic obser- 
vatory is tlie dt?pendabi 1 ity . The electronic auci otlier components used 
in an automatic observatory are inevitably more compK--'. tlian the instru- 
mentation in the conventional system, and hence tliere is mori* chance of 
loss of data arising from various equipmental ma i func t ior s . Let us not 
forget , liowever, how often magnetograms especially from liigh latitude 
observatories are useless or nearly so during disturbances when they are 
needed most. Also, manual operations that are assumed to be dependable 
do not appear to be as reliable as are assumed theoretically. Examples 
proving this point are numerous. Magnetogram traces are often lost 


Ijccause of po(U* optica; woik Lii the record In poor phott>-pro<. ess intd; , 
or careless phutoprapli ic and processing work Cor tlie repnuluct ion of 
the inagne togrnins . f)ther »!xainples of causes for loss of data are; a 
lack of, or errors in, the indent if icat Lons of the tiiree ctmiponents ; 
errors, or uncertainties, regarding the signs for the ordinates in tlie 
inagne tograms ; or, in the case of tabulated hourly values, misreading the 
inagnetogram traces because of reflected images from off-scale excursions. 

In the foregoing discussions I have expressed usefulness of digi- 
tized geomagnetic data and presented arguments advocating automatic 
observatories that will produce digitized data directly from Llie magneto- 
meters. However, wiieii such digital data become available in a machine 
readable form, there may arise a serious problem that lias not been 
experienced previously. This problem concerns the enormoviS advantages 
given to scientists who have direct access to computers. Scientists in 
those areas of tlv' world wliero accessibility to computers is poor may 
in practice be denied opportunities to participate meaningfully in the 
analysis «.>f geomagne t tc data. In those countries in whicii use oi emn- 
puters is limited, geomagnetic observations may become, in essence, a 
service to the scientists in more technologically advanced countries, 
riierehy active interests and sciiuiUtfic participations of the conntritis 
of the former category may gradually decline. This is not a desirable 
circumstance and certainly not in line with the spirit of the lAGA. 

Such 3 trend is coutiurv to the traditions of the international cooper- 
ation Lp the Polar Years, IGY, ICC, IQSY, lASY, and in other international 
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t’lent of an 1 iitor iiat Lona i Co'Uct f»)r (lootna^no L i c Stiidic>. llu- fundai iciii 1 
purposes of tlu Ceiitur are: (i) to nial'U a coinputk. r ‘-y-.t >.Tii , p'otLi.r"'. 

anil other i-qiiipinjnL noeded ii)r data analysis availalili ti> e i eiu i st.-. , 
and (li) to provide oi<i)ortun i t ies Li> Liiose scientists who at theit lu'ina 
institutions havi‘ only limiteil access to such inoden. f a<. i 1 i t i i'.-> . rlu- 
Center should havt' a permancuit staff «^f liiyh qua 1 i I ii .ti iim .o lliai vi.il- 
inj; scientists can riceive advice from them either .sc ieui i i i ea I I' i r 
I <- clini ca 1 ly . Positions of "visiting; professiu's" may provide Si.I».ntiliw 
stii:iiilus. rite Center v.’ould be a lot^ical for a Di; 4 jLal h.iia C».'nL> r 

such as the one tliat has been sn^^ested hy Dr. .1. . Cain. 

1 he liost country for tlie Ciu'.tiT would sei.‘m Li) Dave i.o U'tisly tlu- 
folUwin^ conditions: (i) Tlie computer Lecliru' I o>;y is hi.^ihly aiiv.uued, 

(ii) rhere e.xists -i::able group t/f scientists a». lively onpaved in uieo- 
magiu'tic researclies, (iii) N'o d iscr iminat ion ; made ri'gardin/, tiu 
visitors' nat iona 1 i t i es . 

.Wed less to say, tlie problem of funds lor t!ie estal-' I i shment iHil 
operation i>f such a Center woulil he enormously difficult. However, iii- 
thusiasm and an accumulation of devote’d efforts can often solve' sei'mingly 
insur'nount.ib le tasks in many areas of int-Tnat iona I i-nte rpr ise' . Hut the 
prospect of materialising tlu- Center wituld be sero if ine enthu.'. i a m of 
the scientists with wltitm tlie matter is concerned is onl" lukewarm. I 
believe tiiat tlie key question is wliether sucli a (.’enter i.s needc d or merely 
de .i i rab le . if the latter, I see no prospect of success and any effort 
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will hc! in vain. Wlio then is to answer this question? Tlie an.swer must 
of course come from all concerned; hut those who would benefit most from 
the plan must provide the driving force. It v^^ou Id be an idle thought if 
one supposes that a country, or countries, in which modern facilities 
art abound would initiate such a drive. The prospect of success will 
ttnis depend first on the strength of the motivation of those who need 
the Center and secondly on tiie understanding of their need on the part 
of those who can assist financially. 
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